The spin and charge gaps of the half-filled N-leg Kondo ladders 
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In this work, we study N-leg Kondo ladders at half-filling through the density matrix renormal- 
ization group. We found non-zero spin and charge gaps for any finite number of legs and Kondo 
coupling J > 0. We also show evidence of the existence of a quantum critical point in the two 
dimensional Kondo lattice model, in agreement with previous works. Based on the binding energy 
of two holes, we did not find evidence of superconductivity in the 2D Kondo lattice model close to 
half-filling. 
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Recently, a great deal of excitement has been gener- 
ated with the discovery of several heavy fermion com- 
pounds which exhibit a quantum critical point (QCP)1*£. 
Several anomalous behaviors have been observed close 
to the QCP and one of the main problems of this area 
is to try to understand this puzzle. The Kondo lattice 
model (KLM) is the simplest model believed to describe 
heavy fermion materials^,. This model incorporates the 
interaction between the localized spins and the conduc- 
tion electrons via an exchange interaction of strength J. 
This interaction is believed to be the most important 
in heavy fermion compounds at low temperatures. At 
small values of J, due to the indirect Ruderman-Kittel- 
Kasuya-Yosida (RKKY) interaction, antiferromagnetic 
long-range order (LRO) is expected, whereas for large 
J a paramagnetic state emerges. The competition be- 
tween these two states leads to the QCP, as first pointed 
out by DoniacbA Many theoretical works have concen- 
trated on the one-dimensional (ID) model, which can be 
treated more easily. In particular, at half-filling, several 
numerical calculations show that the ID KLM is always 
a gapped spin liquid with no magnetic quasi-LRO MiT . 
Away from half-filling, there is power law decay of spin 
correlations, for small values of J (see, however, Ref. @). 
On the other hand, for the two-dimensional (2D) KLM 
athalf-filling, several studies support the existence of a 
QCPAiLii (see also Ref. 

In this work, we present the spin and charge gaps of the 
N-leg Kondo ladders at half-filling calculated by the den- 
sity matrix renormalization group (DMRG) technique^. 
The N-leg Kondo ladders consist of N Kondo chains of 
length L coupled by the hopping term. The 2D system 
is obtained by taking both N and L to infinity. As far 
as we know this is the first study of the iV-leg Kondo 
ladders. We will also show evidence of the existence of 
a QCP in the 2D Kondo model at half-filling in agree- 
ment with previous studies. Moreover, our results for the 
binding energy of two holes do not show any evidence for 
an effective hole-hole attraction in the 2D KLM for any 
Kondo coupling. All the results we will report here cor- 
respond to the conduction electron density n = 1. It is 
important to note that for this density, techniques other 
than the DMRG can be applie dSi 1 ^ 11 . However, away 
from half-filling, methods like Quantum Monte Carle 11 



are plagued by the "sign problem". Therefore, our re- 
sults at half-filling should serve as benchmarks studies 
before get results at more general fillings. DMRG re- 
sults away from half-filling are in progress and will be 
presented elsewhere!^. 

We considered the N-leg Kondo Hamiltonian at half- 
filling with M — NxL sites 

<i,j>,<r j 

where Cj a annihilates a conduction electron in site j 
with spin projection a, Sj is a localized spin \ opera- 
tor, Sj — | 2 a/3 c ] a^apCj p is the conduction electron 
spin density operator and ar a p are Pauli matrices. Here 
< ij > denote near-neighbor sites, J > is the Kondo 
coupling constant between the conduction electrons and 
the local moments and the hopping amplitude was set to 
unity to fix the energy scale. We investigated the model 
with the DMRG technique with open boundary condi- 
tions. We used the finite-size algorithm for sizes up to 
N x L = 40, keeping up to m — 2000 states per block. 
The discarded weight was typically about 10~ 4 — 10 -8 in 
the final sweep. Although the DMRG is based on a one- 
dimensional algorithm, it has been applied to low dimen- 
sional systems. The procedure consists in mapping the 
low-dimensional model on a ID model with long range 
interactional^^. Several important questions, such as 
the existence of striped phases and/or superconductivity 
in the 2D t — J model, have been addressed with this 
procedure!!. Here, we used this same method to study 
the N-leg Kondo ladders. 

Our first concern is to check the applicability of the 
technique to the N-leg Kondo ladders. We can get some 
insight by looking at the ground state energies. We have 
checked that, for small slightly doped clusters with up 
to 12 sites, we can reproduce with the DMRG the same 
energies obtained independently by exact diagonalization 
(ED). As we increase the system size the Hilbert space 
increases exponentially and ED is no longer a suitable 
method to treat the system. This is the main advantage 
of the DMRG. In order to illustrate this advantage, we 
will show that we can obtain the ground state energies 
of the N-leg Kondo ladders in a controlled waj^. 
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2x10 


4x4 


m 


J=0.0 


J=2.0 


J=0.0 


.7=0.3 


J = 3.0 


100 


-27.56929 


-39.15411 


-21.85340 


-22.12936 


-41.11699 


200 


-27.62807 


-39.15909 


-21.88557 


-22.19189 


-41.25039 


300 


-27.63982 


-39.15988 


-21.88816 


-22.25239 


-41.28139 


400 


-27.64290 


-39.16007 


-21.88845 


-22.29403 


-41.29605 


600 


-27.64407 


-39.16015 


-21.88853 


-22.33682 


-41.31698 


800 








-22.36031 


-41.32227 


1200 








-22.38726 


-41.32606 


oo 


-27.6455 


-39.1602 


-21.8885 


-22.4413 


-41.3401 


exact 


-27.64486 




-21.88854 







Table I: The ground state energies of the Kondo laticce model 
for the clusters 2x10 and 4x4 as a function of the number of 
states m kept in the DMRG truncation for some Kondo cou- 
plings. The energies at m = oo were obtained by an extrap- 
olation, see text. 



In table 1 we present the ground state energies as a 
function of the truncation m for the 2-leg (4-leg) Kondo 
ladders with length L = 10 (L = 4), The energies for 
m = oo were obtained by an extrapolation assuming a 
dependence of the type e(m) = eoo + a/m + b/m 2 . As 
we can see, for the J = (free case) with a small num- 
ber m of states the relative error is smaller than ~ 10~ 4 . 
This free case result does not appear to depend on the 
number of legs, as shown for the 2- and 4-leg ladders. 
However, a small coupling J makes the wave function 
much more complex and a larger number of states is 
necessary to keep the same accuracy. Furthermore, in 
this weak coupling regime, we need a much larger num- 
ber of states (and consequently a longer computational 
time) to keep the algorithm stable as the number of legs 
is increased. An increase in computer effort, probably 
exponential, such as this has also been observed in other 
system9i*!*A&. This was the main reason why we restricted 
our calculations to N < 6 and large values of J. On the 
other hand, in the strong coupling limit, a smaller num- 
ber of states is necessary to keep the same accuracy. For 
J>0 the conduction electrons form unbreakable mobile 
singlets with the localized spins, effectively behaving like 
free spinless fermions. 

Having set these benchmarks, we are now able to cal- 
culate the spin and charge gaps of the N-leg Kondo lad- 
ders with M = NxL sites at half- filling. These are 
defined as A S (L) = e(M,l) - e(M,0) and A C (L) = 
e{M — 2,0) — e(0,0), respectively, where e(N c ,s) is the 
ground state energy of the system with N c conduction 
electrons in the subspace of total spin si^L. In Fig. QJa) 
we present the spin and charge gaps for the 2-leg Kondo 
ladder with length L = 10 (open symbols) and L = 20 
(closed symbols) as a function of the Kondo coupling J. 
As we see, finite size effects are very small. In the strong- 
coupling limit the spin and charge gaps increase linearly 
with J, and can be estimated by second order perturba- 
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Figure 1: (a) The spin and charge gaps for the 2-leg Kondo 
ladder with L = 10 (open symbols) and L = 20 (closed sym- 
bols). The circles (squares) correspond to the spin (charge) 
gap. Inset: zoom of the region with < J < 2. (b) Same as 
(a), but for the 3- leg Kondo ladder with L = 8 (open symbols) 
and L — 14 (closed symbols). 

tion theory with respect to i/J— . At this order, the spin 
and charge gaps of the A-leg Kondo ladders are given, 
respectively, by A s 2) = J - and A^ 2) = 3 J/2 - 

27^ + 7^(^-1)^, where j N = l-coa(7rJV/(JV+l)). 
We used the infinite-size algorithmic to determine the 
spin and charge gaps of the 2-leg Kondo ladders in the 
strong coupling limifeSS. Our numerical calculations are 
in very good agreement with these results. For example, 
for J = 40 we obtained A s = 39.752 and A c = 57.057, 

(2) 

very close to the perturbative results As = 39.750 and 

(2) 

A y c = 57.037. On the other hand, for small values of J 
the spin gap has a different dependence on the coupling 
constant. 

The results we have found for the 2-leg Kondo ladder 
are very similar to those of the ID KLMt In this case, 
for large J the spin gap goes linearly with J, while for 
small J it decays exponentially as A s ~ exp(—l/anJ), 
where a ~ 1.41&. This kind of behavior does not seem 
to depend on the number of legs. The spin and charge 
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gaps of the 3-leg Kondo ladder show similar behavior, as 
shown in Fig. db). Based on the spin and charge gap 
behaviors of the 1-, 2- and 3-leg Kondo ladders, we expect 
non-zero spin and charge gaps for the half-filled TV-leg 
Kondo ladders. It is interesting to note that the spin gap 
has similar behavior whether the number of legs is even 
or odd. The same does not happen in TV-leg Heisenberg 
ladders, where the spin gap is zero for odd TV and different 
from zero for even TV— (for a review see Ref. l22T l. This 
difference arises from the fact that the Kondo Ladders 
are coupled just with the hopping term, differently from 
the Heisenberg Ladders. So, the same analysis that is 
done in the N-leg Heisenberg ladders (see Ref. l22Tl . does 
not apply in the N-leg Kondo ladders. Moreover, we 
might naively expect that the inclusion of the Coulomb 
interaction U in the N-leg Kondo ladders would lead to 
a spin gap with a distinct dependence on the number of 
legs (at least for U ^S> ), as in the Heisenberg ladders. 
However, most probably, even the presence of U will not 
lead to this distinct behavior. Note that the ID Kondo- 
Hubbard model at half-filling has a non-zero spin gapl&, 
differently from the ID Heisenberg model. 

We would still like to determine whether the spin 
gap remains finite or not when TV goes to infinity (af- 
ter L — ► oo). Approximate approaches&ifi and Quantum 
Monte Carlo results^ support the existence of a QCP at 
J c /t ~ 1.45 ± 0.05. Below this value antiferromagnetic 
LRO appears with a zero spin gap. It is a very tough 
task to determine the presence of LRO directly by the 
value of the spin gap, since for small values of J the gap 
is very small, and as we increase the number of legs we 
lose accuracy in the energies. 

Instead of looking at the spin gap to check for LRO, we 
can calculate the Fourier transform of the spin-spin cor- 
relation function S (q) = jgEf 1 , r - 2 e ' ' ri ~ r2 ' (Sf*i ' Sf 2 ). 
In Fig. EI a), we show the spin structure factor S (q) for 
the 2-leg Kondo ladder as a function of momentum q x 
for J = 0.5. The peak at q = (tt,tt) is a signature 
of some kind of antiferromagnetic order, not necessarily 
long-ranged. For small values of J we expect this peak to 
increase as we decrease J, since the RKKY interaction 
starts to dominate over Kondo compensation. Indeed, 
this is what we have seen. In Fig. EJb) , we present the 
spin structure factor S (q — (tt,tv)), for the clusters 2x10 
(squares) and 3x8 (circles), as a function of the Kondo 
coupling J. Note the divergence as J — * 0. This kind 
of behavior was also observed in the magnetic suscepti- 
bility of the one-dimensional KLM^. Antiferromagnetic 
LRO exists if liniM^oo S (q= (ir, tt)) /M is non-zero. In 
Fig. E{c), we present S(q= (it, tt)) /M as a function of 
VM for several values of the Kondo coupling. As can be 
seen, even working with small clusters, our results sug- 
gest that for J > 1.6 S (q = (tt, it)) /M tends to zero as 
we increase M, while for J < 1.3 it tends to a finite value. 
This result is compatible with J c — 1.45 as suggested by 
other method aW 1 . 

Finally, we would like to check for evidence of super- 
conductivity in the 2D KLM close to the QCP. Several 
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Figure 2: The spin structure factor S (q) for the N-leg Kondo 
ladders, (a) S (q) vs momentum q x for the 2-leg Kondo lad- 
der with length L — 10 and J = 0.5. (b) S (q— {it, 7r)) vs 
the Kondo coupling J for the clusters 2x10 (squares) and 3x8 
(circles), (c) S (q= (n, tt)) jM vs 1/M for several Kondo cou- 
plings. The sizes of the clusters used were M =2x2, 3x4, 4x4, 
4x6, and 6x6. 




Figure 3: The binding energy Ab for the 1-, 2- and 3-leg 
Kondo ladders as a function of the Kondo coupling J. 



heavy fermions compounds, like Celn 3 and CePd2Si2, can 
be tuned to an antiferromagnetic QCP by application of 
external pressure^. Superconductivity was observed at 
low temperatures close to the QCP of these materials. 
Our aim will not be to find proof of superconductivity at 
finite hole doping close to a quantum critical point in the 
2D KLM, which would entail the calculation of the pair 
correlation function. Instead, we can gain some insight 
by looking at the binding energy of two doped holes, de- 
fined as A B = E{2) + £7(0) - 2£7(1), where £7(n) is the 
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ground state energy with n holes. If the holes are bound 
Ab < 0, and this is indicative that attractive effective 
forces are present. This is a necessary, though not suf- 
ficient, condition for a possible condensation of the holes. 
In Fig. El we present the hole binding energy for the 1-, 2- 
and 3-leg Kondo ladders. In the thermodynamic limit the 
binding energies vanish if holes do not form bound states 
and is asymptotically positive^. As we see from this fig- 
ure, this is exactly what we have observed. Increasing the 
cluster from 2x10 to 2x20 the binding energies decrease. 
As we see, the holes do not bind for any Kondo coupling. 
Similar results were also observed in the Kondo chairA 
These results suggest that, at least close to half-filling, su- 
perconductivity will not emerge in the 2D Kondo lattice 
model. The fact that we have not found bind holes also 



suggest the the KLM may not be the minimum model to 
describe correctly the Heavy Fermion compounds. 

In summary, we presented, for the first time, a system- 
atic study of the 2- and 3-leg Kondo ladders at half-filling 
by using the DMRG. Based on the behavior of the spin 
gap of the 1- to 3-leg Kondo ladders, we expect non-zero 
charge and spin gaps in the N-leg Kondo ladders for any 
J > 0. Using small clusters we were able to determine 
the presence of a QCP, in agreement with previous works. 
Our results also suggest that there is no superconductiv- 
ity close to half-filling in the 2D Kondo lattice model. 

The author thanks E. Miranda, E. Dagotto, P. G. 
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